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Abstract 

We derive, via fourth order perturbation theory, an expression for the Joseph- 
son current through a gated interacting quantum dot. We analyze our expres- 
sion for two different models of the superconductor-dot-superconductor (SDS) 
system. When the matrix elements connecting dot and leads are featureless 
constants, we compute the Josephson coupling Jc as a function of the gate 
voltage and Coulomb interaction. In the limit of a diffusive dot, we compute 
the probability distribution P{Jc) of Josephson couplings. In both cases, vr 
junction behavior (Jc < 0) is possible, and is not simply dependent on the 
parity of the dot occupancy. 

Magnetic impurities in Josephson junctions tend to degrade the critical current This 
result has been known since the work of Kulik |]T|, who analyzed a simple model in which 
spin-preserving as well as spin-flip hopping processes (with amplitudes t and t', respectively) 
are present. The critical current in such a junction is given by the expression Jc = J^^ • 
(|tp — |t'P)/(|tp + It'l"^), where J^^ is the Ambegaokar-Baratoff critical current. Since spin- 
flip tunneling results in a sign change of the Cooper pair singlet (t| to |t), the spin-flip 
hopping contribution is negative, and reduces Jc |]5[. When Jc < 0, one has a vr junction, for 
which the ground state energy is minimized when the superconducting phase difference is 
(5 = TT. A ring containing a single vr junction will enclose trapped flux 0. A related effect 
occurs in Josephson tunneling through a ferromagnetic layer ||8|-[ll|, and recent experiments 
on Nb-Cu^Nii_a;-Nb junctions suggest that Jc < states have been observed 



In this paper, we investigate Josephson coupling mediated by a quantum dot [T^, gen- 
eralizing the case of a single impurity to a system with many quantized energy levels. We 
derive flrst a general expression, within fourth-order perturbation theory, for the Josephson 
coupling Jc (Jc = 2eJc/h). We then consider two models for the tunneling amplitudes taj 
from the leads to the dot. As we shall see, in contrast to the single impurity case, the 
parity of the number of electrons on the dot, A^o; does not uniquely determine the sign of 
the Josephson coupling. 

It is well-known in the theory of elastic co-tunneling that if the Coulomb repulsion U 
is large then the conductance through a dot is independent of U. Correspondingly, we flnd 
the critical current is insensitive to U in this regime, and furthermore when the system is 
close to a charge degeneracy point, the probabihty distribution P(Jc) for a disordered dot 
has universal properties. 

Our Hamiltonian is a sum of three terms: 
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Here, Sct,d arc the areas of the a electrode and the dot, V is the gate voltage on the dot 
and U is the Hubbard interaction on the dot. Electron states on the dot are assumed to be 
disordered by a spin-independent random potential. 

We calculate the Josephson current via fourth order perturbation theory in T^tun- For 
this to be applicable both the typical distance between consecutive energy levels Se and the 
gap in the superconductors must be large compared to the broadening of the individual level 
due to tunneling: <C Ss,Aa, where Fq = 7ii'a{\ta\'^) / ^d, where is the metallic density 
of states per unit area in the a electrode. 

We also have to be far enough from charge degeneracy points, where the gap for charge 
excitations on the dot vanishes. We measure V relative to the charge degeneracy point, 
which is equivalent to setting fi = sij^^^^i + UNq. 

In computing the fourth order correction to the ground state energy, we only consider 
terms which depend on the phase difference 5 between the two superconductors. 



Jee = [( Ai cosh ei + ^ + V){C + C' + 2V + U) ( A2 cosh 62 + ^' + V)]'^ 

+ [( Ai cosh Oi + ^ + V) ( Ai cosh 61 + A2 cosh 62) ( A2 cosh 62 + ^' + V)]'^ 
+ [(Ai cosh ei + C + V){C + C' + '^V + U){A2 cosh 02 + C + V)r^ 

Jhh = [(Ai cosh 0i-C-V + U){-C-C' -'^V + 3C/)(A2 cosh e2-C'-V + U)]'^ 

+ [( Ai cosh Oi-^-V + U)(Ai cosh 61 + A2 cosh 62) ( A2 cosh 62 - ^' - V + U)]'^ 
+ [( Ai cosh Oi-^-V + U)(-^-^' -2V + 3U) ( A2 cosh e2-^-V + U)r^ 

Jeh = [( Ai cosh 0i-C'-V + U) ( Ai cosh 01 + A2 cosh 02 + C- ( A2 cosh 02-C' -V + U)]'^ 
+ [( Ai cosh 01 + ^ + V) ( Ai cosh 01 + Ai cosh 02 + ^- C') ( A2 cosh 02 + ^ + V)r^ 
+ [( Ai cosh 01 + ^ + V) ( Ai cosh 0i + A2 cosh ^2) ( Ai cosh 0i-^'-V + U)]~^ 
+ [( A2 cosh ^2 + { + V^) ( Ai cosh 0^ + A2 cosh ^2) ( A2 cosh 02-^' -V + U)] 
+ [( Ai cosh 0i + ^ + V) ( Ai cosh 0^ + A2 cosh 02 + ^- O ( Ai cosh 0i-^' -V + U)] 
+ [( A2 cosh 02 + C + V) ( Ai cosh 01 + A2 cosh 02 + C- ( A2 cosh 02 - C' - V + U)]-' , (3) 



Ej{S) = -4i/ii/2AiA2Ree^^ ^ ti,ti/t;,%^ee(e„ C/) 





(2) 



where Jcc, Jhh and Jch ^.tg given by 
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where ^ = e — fi. 

To proceed further we have to make some assumptions about the tunnehng amphtudes 
and energy spectrum of the dot. Initially, we ignore any structure to the hopping amplitudes 
and set t^j = ta- This choice of the amplitudes is equivalent to a situation where both leads 
are connected to the same point on the surface of the dot. 

We further assume the energy spacings on the dot satisfy 5e ^ Aq,. The summation 
over states on the dot can then be recast as an integral. We obtain Ej = Jc{l — cos 5), with 

Jc = ^i^l^i^2^i^2 dx dy Jcc{B^x,Bey) 

TT^ JO Jo L 

+Bl Jhh(-5ha;, -B^y) - B^B^, J^^{B^x, -B^^y)] , (4) 

where -Be,h are the distance to the top/bottom of the dot spectrum from the Fermi level on 
the dot. 

By numerical integration we obtain JdV)-, shown in fig. la for different values of U . For 
small values of U the Josephson coupling always is positive. However, above some critical 
value Uc there is a finite interval of voltages in which Jc < 0. The phase diagram of the 
junction is presented in fig. Ib-c. Note that our approach does not apply close to the V = Q 
otV = U lines where the Coulomb blockade is lifted and the dot ground state is degenerate. 

Disordered dot - When the leads are connected to the dot at two different points, it is 
necessary to account for the effects of disorder. In this case, we have taj = ta '?/'j(Ra), where 
'?/'j(r) is the wave function corresponding to energy level Sj on the dot. These wave functions 
are functionals of the disorder potential. The Josephson coupling is now a random quantity 
and we must find its distribution, which is possible since the statistical properties of the tpj 
are well-studied. 

Let us define the dimensionless local density of states, 

^-D T3 ^ <^<^(R'l>R-2) - G^(R.l,R-2) 

Pi^l-tll, -tl2j — — 

ZTTZZ/d 

= z/,i^^,(Ri)^;(R2)5(^-^,) , (5) 
j 

where v is the dot density of states u^ = Ij {SaSs), and G^'^ are the usual retarded/ advanced 
single particle Green's functions. 



As it was discussed in [14| for u less then Thouless energy E'xh these objects have very 



simple disorder-averaged statistical properties: 



all pairwise 
contractions 



implying that the function is distributed according to the functional 

PHocexpj- rdu-^] . (7) 



The Josephson energy is then a bilinear functional of 
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+ / / dujduj'p^p^iJi,^{uj,uj') - / dujduj'p^p^:J^^{uj,uj')) (8) 



— oo J — oo 
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where we define dimensionless conductances Qa = i^ai^dl^a] ■ 

Formally, our task is straightforward. We have to obtain critical current distribution 
function P(Jc) = {^{Jc — Jc[Pui\))- We do that as follows. First, we find the eigenvalues of 
the matrix 7(^7, cu'), where 

{Jee(co', Uj') + Jeeiyj' , ^) if Co" > 0, tu' > 
-Jeh(LL'',u;) if w < 0, cu' > 

Jhh(co', t^') + Jhh{'^\ ^) if ci^ < 0, uj' < 0. 

The distribution P{Jc) now has the form 

P{Jc) = {6 (Jc-4AV^2 [p-J-p]] 

du _,-„j„-|-|- .... s .9 ,\~i/2 

A 
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/ ^e"^"'^ n (l - 4m(fe) A^^i^^ a) , (10) 



where the product is over all eigenvalues A of the matrix J. Normalization determines the 
appropriate branch cut. 

If both leads are connected to the same point on the surface of the dot, then oc 
Jdu Jduj' J{u,uj'), which we have previously shown to be negative for certain V,U. Hence 
J is not in general positive-definite, and we anticipate some finite probability for tt junction 
behavior. 

Our numerical procedure consists of several steps. We introduce a cut-off W and confine 
UJ to the interval —W/2 < oj,uj' < W/2. Further, J{u!,u!') is put on a uniform frequency 
grid. The Josephson coupling now is Jc = 4,/S?gig2 Y.s,s' Jss'PsPs' where the 2N x 2N matrix 
Jss' = {W/2N) J{ijJsUj's) and ps = {W/2NY^'^ p^/, {PsPs') = 5e 5ss'- Jss' is diagonalized and its 
eigenvalues A^ are substituted into ([Tol). To make this integral more suitable for numerical 
evaluation a contour transformation is performed: 

nJc)=j^ — e^"'^^n(|lT4«fcA2^i^72A.|) ' sin(-^e(TnA,-l)j. (11) 

One should pick upper sign for Jc > and lower sign for < 0. 

The results of numerical calculation are quite insensitive to W and A^. The distribution 
function P{Jc) is presented on fig. 2. It has a shape of asymmetric bell and it is non- 
zero for Jc < per previous discussion. The total probability of having vr-junction: P^^ = 
P{Jc) dJc, is shown in fig. 3. 

Universal Limit - Generally the distribution P{Jc) of Josephson couplings is dependent 
on nonuniversal parameters such as V and U. In the limit f/ ^ Ai = A2 = A, though, 
there are universal aspects to the distribution. If the dot is close to the charge degeneracy 
point, i.e. r<^y^Aorr^f/ — l^^A then the matrix J is independent of both U 
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and V. Restricting our attention to the first case, we find the matrix J is still given by 
with Jhh = and 

/•oo /"OO 

Jee ^ dOi (162 [(A cosh ei+uj){A cosh 6*1 + A cosh 6*2) (A cosh 62 + 00')] 
Jo Jo 

Jeh ~ / ^611 / 6/6*2 [(A cosh 6*1 + cj) ( A cosh 6*1 + A cosh 6*2 + u; - u') (A cosh 6*2 + cj)] " . 
Jo Jo 

Thus, J{lj,u!') = F{u!/ A,u!' / A)/ where F{x,x') is a universal function of its argu- 
ments. We then find -P(Jc) = Jab fi-^c/ Jab), where f{x) is universal and Jab = TT{G)A/h. 



The average conductance through the dot at V = A was calculated in ||T^ to be (G) = 
{27ieyh)g^g2 (fe/A). 

There are several consequences of this universality. First, it means that all moments of 
this distribution are proportional to powers of Jab, with universal coefficients. In particular, 
the ratio of the RMS and mean critical currents is ((A 7^)^) / ( Jc) ~ 1.59. Second, the 

probability of having a vr-junction is a universal number: P^, = ^l^dxf {%) ~ 0.19. 

Odd Nq - When the parameters of the dot are chosen in such a way that there is a single 
level '0' occupied by only one electron, equation must be modified. Consider a Cooper 
pair tunneling from the left to right superconductor. If none of two electrons passes through 
'0' then their contribution to Ej is already included in (|^). Only those events where one or 
both electrons tunnel through '0' will modify the expression for Ej. We find the corrections 
are given by 

AJee(u;,Cj') = SeSiu') (^Jee{uJ,0) - ^J,^{iJ,0)^ 9^ 

AJhh(a;,cu') = SeSico') (^Jhh(w,0) - ijeh(0,cu)) Oi-co) 

A Jeh(u;, uj') = (fe)' 6{uj) 6{uj') Jhh(0, 0) . (12) 

In the case of a single impurity level, we recover previous results Overall, this electron 
increases our chances to get vr junction. It is clear, however, that the influence of that 
electron is small for metallic dot where 6e <^ A. 

Discussion - Two superconductor-dot-superconductor models have been considered 
here. Assuming that both leads are attached to the same point on the dot surface and 
thus, ignoring any structure to the hopping matrix elements taj, we find a critical U above 
which Jc can be driven negative by an appropriate dot gate voltage, V. As one can see 
from d^) and there are six terms giving positive contribution to Jc and six terms giving 
a negative contribution. They each depend in different ways on the gate voltage V, and one 
can suppress positive terms by varying V. 

Another feature of our result is that it is derived for an even number of electrons on 
the dot. In the case of an Anderson impurity, the sign of Jc depends on the parity of the 
impurity level occupancy Nq [Q. Rather than being a parity effect, we offer the following 
interpretation: If two electrons of the Cooper pair both tunnel through empty states or 
both tunnel through filled states this gives positive contribution to Jc (first and second 
term of (Q)). To obtain a negative contribution (last term of (0)) one electron has to pass 
through filled states and the other through empty states. Thus, for a completely filled 
(empty) Anderson impurity both electrons must pass through filled (empty) states, yielding 
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a positive Jc. For a singly occupied orbital, one electron always tnnncls through empty 
state while the other tunnels through filled state, and J^. is negative. Therefore, the sign 
and magnitude of Jc is determined in part by phase space considerations. 

In the second model we allow leads to be 'connected' at different points on the surface 
of the dot. As a consequence the disorder has to be treated properly. We calculated the 
distribution function of critical current for the ensemble of junctions for different parameter 
values. The probability of tt junction as calculated from this distribution has very reasonable 
values to expect that such tt junction can be found experimentally. It is also predicted 
that the distribution function should possesses a remarkable universal property in the limit 
[/> A. 

There are questions that remain unanswered. Within the framework of our first model, 
in order to have tt junction the value of U must exceed some critical value Uc- From Fig. 3 
one might infer that ^4 > for the second model as well. It has to noted, however, that our 
numerical procedure does not rule out the possibility of exponentially small tail for [/ — > 0. 
If this is true then corresponds to a crossover rather then a phase transition. 
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FIGURES 





FIG. 1. (a) JciV) versus V/U for U = 5, = 80, Bi, = 20 (solid); U = 2, B^ = 33.3, 
-Bh = 66.7 (dashed); U = 3, Bg = 33.3, B^ = 66.7 (dot-dashed), (b) phase diagram for B^. = 50 
(soUd), Be = 33.3 (dashed), and B^ = 10 (dot-dashed) with B^ = 100 — B^ The tt phase hes above 
the curve in each case. 




FIG. 2. Probabihty distribution P(Jc) ior U = 3, V = 0, W = AO (soUd) and U = 3, V = 1.5, 
and ly = 40 (dashed). = 100 in both cases. The area under the total curve (i.e. out to | Jd = cx) 
is unity in both cases. 
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FIG. 3. versus U/A for V = 0.1 A. For very large values of C//A, 0.19. 
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